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ABSTRACT 

We report on the Ka-band (26—40 GHz) emission properties for 10 star-forming regions in the nearby 
galaxy NGC 6946. From a radio spectral decomposition, we find that the 33 GHz flux densities are 
typically dominated by thermal (free-free) radiation. However, we also detect excess Ka-band emission 
for an outer-disk star-forming region relative to what is expected given existing radio, submillimeter, 
and infrared data. Among the 10 targeted regions, measurable excess emission at 33 GHz is detected 
for half of them, but in only one region is the excess found to be statistically significant (« 7a). We 
interpret this as the first likely detection of so called 'anomalous' dust emission outside of the Milky 
Way. We find that models explaining this feature as the result of dipole emission from rapidly rotating 
ultrasmall grains are able to reproduce the observations for reasonable interstellar medium conditions. 
While these results suggest that the use of Ka-band data as a measure of star formation activity in 
external galaxies may be complicated by the presence of anomalous dust, it is unclear how significant 
a factor this will be for globally integrated measurements as the excess emission accounts for <10% 
of the total Ka-band flux density from all 10 regions. 

Subject headings: radio continuum: general - dust - galaxies: individual (NGC 6946) 



1. INTRODUCTION 

Accurate separation of Galactic foreground emission 
in cosmic microwave background (CMB) experiments 
remains a major challenge in observational cosmology. 
Besides the standard foreground components, free-free, 
synchrotron, and thermal dust emission, an unknown 
component has been found to dominate over these at 
microwave frequencies between ~ 10 — 90 GHz, and 
is seemingly correlated with 100 p,m thermal emis- 
sion from interstellar dust (e.g. iKoeut et al.~l 119961 : 
Ide Oliveira-Costa et al. Ill997t ILeitch et al. 1119971 ). 

The first detection of this 'anomalous' dust-correlated 
emission was attributed to free-free emissio n from shock- 
heate d gas in sup ernova (SN) remnan t s (ILeitch et al. I 
Il997h . However, iDraine &; Lazariaril (|1998af ) demon- 
strated that the observed microwave excess could not 
be due to free-free emission from hot gas, as this sce- 
nario would require an unrealistically high energy injec- 
tion rate (i.e. >2 orders of magnitude larger than that 
from SNe), and suggested two alternative models; electric 
dipole rotational emission from ult rasmall (a < 10~ 6 cm) 
grains (|Draine fc Lazarian Ill998bf ). and magnetic dipole 
emission from thermal fluc tuations in the magnetiza tion 
of interstellar dust grains (jDraine fc Lazarian II1999T ). 

At present, the most widely accepted explanation 
for the anomalous dust emission is the spinning dust 
model in which rapidly rotating very small grains, 
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having a non-zero electric dipole moment, produce 
the ob served microwave emission (jDraine fc Lazarian I 
1998b] : This model was r e cently refined by 



Ali-Hai'mo ud, Hirata. k. Dickinson"! (|2009f ) who imple- 



mented the Fokkcr-Planck equation to compute grain 
angular velocity distributions using updated grain op- 
tical properties and size distributions. To date, anoma- 
lous dust emission has not been directly detected out- 
side of the Galaxy. In this Letter we present the first 
likely detection for an outer-disk star-forming region in 
the nearby galaxy NGC 6946. 

2. OBSERVATIONS AND DATA ANALYSIS 

We targeted 10 star-forming regions in the 

nearby galaxy NGC 6946 (d » 6.8 Mpc; 

iKarachentsev. Sharina. fc Huchtmeier I [2000) . including 
its starbursting nucleus, chosen due to their wealth of 
existing and forthcoming mid- and far-infrared spec- 
troscopic data collected as part of t he Spitzer Infrared 
Near by Galaxies Survey (SINGS; iKennicutt et. al. 1 
2003) and the project Key Insights on Nearby Galaxies: 
a Far-Infrared Survey with Herschel (KINGFISH; PI. 
R. Kennicutt). The location and identification number 
of each region is overlaid on the stellar continuum 
subtracted 8 /im image, which traces emission from 
very small grains (VSGs) and polycyclic aromatic 
hydrocarbons (PAHs) , in the top panel of Figure [T] 

2.1. Radio Data 

Observations in the Ka-band (26 — 40 GHz) were taken 
on 2009 March 21 using the Caltech Continuum Backend 
(CCB) on the 100 m Robert C. Byrd Green Bank Tele- 
scope (GBT). The CCB simultaneously measures the en- 
tire Ka bandwidth over 4 equally spaced frequency chan- 
nels and synchronously reads out and demodulates the 
beamswitched signal to remove atmospheric fluctuation 
and/or gain variations. Reference beams are measured 
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by nodding 1 f 3 away from the source, and their positions 
are identified on the 8.5 GHz radio map in the bottom 
panel of Figureffl The average FWHM of the GBT beam 
in the Ka-band was «25"among our sets of observations, 
and is indicated by the size of each circle in Figure [TJ 
at the distance of NGC 6946 this projects to a physical 
scale of « 0.8 kpc. 

A detailed description on the performance of the 
CCB receiver, along with information on the reduc- 
tion pipeline a nd e rror estimates, can be found in 
IMason et. al~l ((2009). In addition to any systematic 
errors, a calibration error of 10% was assigned to the 
flux density measur ement at each frequency channel 
(|Mason et. al. II2009T) . To analyze the full radio spectral 
energy distributions (SEDs) of these regions, we com- 
piled multifrequency radio data from the literature. The 

1.4 GHz radio map (14" x 12'.'5 beam) comes from the 
Wester bork Synthesis Radi o Telescope (WSRT)-SINGS 
survey (jBraun et al. I 120071 ) while the 1.5, 1.7, 4. 9, and 

8.5 G Hz data (15" x 15" beam) all come from iBeck I 
(|2007ft . The 4.9 and 8.5 GHz radio data included single- 
dish measurements. 

2.2. Infrared and Radio Photometry 

All infrared data used here were included in the SINGS 
fifth data release. We make use of 8, 24, and 70 /im 
Spit zer data, as well as 450 and 850 /zm SCUBA data 
(see lDale et al. II2005D . to compare the infrared and ra- 
dio properties of each region as these data are at higher 
spatial resolution than the GBT data. Regions 2, 3, and 
5 were not covered by the SCUBA 450 or 850 /im maps, 
while regions 4 and 8 were near the edges of the maps. 
The 8 /im map was correcte d for stellar l i ght u sing a 
scaled 3.6 /an image following iHelou et al. I ipOOl . 

Photometry was carried out on all radio and infrared 
maps after first cropping each image to a common ficld- 
of-view and regridding to a common pixel scale. To accu- 
rately match the photometry of our images to the GBT 
measurements, maps were convolved to the resolution of 
the GBT beam in the Ka-band. The Spitzer data were 
first resolution-matched to the 70 /xm beam using custom 
kernels provided by the SINGS team. Since the Spitzer 
PSFs suffer from significant power in their side lobes, 
we convert to Gaussian PSFs using a standard CLEAN 
algorithm; the Spitzer PSF was removed and the image 
was restored with a Gaussian having a FWHM of 25". 
Flux densities at each wavelength were then measured 
by taking the surface brightness at the location of each 
GBT pointing and multiplying by the effective area of 
the beam. The results from our radio photometry, along 
with 1 a error bars, are plotted in Figure [2] for all regions 
except extranuclear region 4, which are shown in Figure 
® 

2.3. Additional Considerations for the GBT Data 

Given that the beam size varies over the full Ka-band, 
we applied a correction to the flux densities at each fre- 
quency channel as if their beam were 25". Scaling fac- 
tors were derived by computing the photometry on the 
8.5 GHz map using the beam sizes for each frequency 
channel averaged over the entire run. These correction 
factors are 0.89, 0.95, 1.00, and 1.04 at 27.75, 31.25, 
34.75, and 38.25 GHz, respectively. The averaged Ka- 
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Fig. 1. — Top: The stellar continuum subtracted 8 fim map of 
NGC 6946 with the ID's of each targeted st ar-form i ng; reg ion. Bot- 
tom: The 8.5 GHz (15"beam) image from IBeckl 1 120071) showing 
the location of the reference beam positions for each observation. 



band flux density, weighted by the errors from each chan- 
nel over the full band, are given in Table 1 along with 
uncertainties; the corresponding effective frequency is 
«33 GHz for each position (i.e. 32.83 ± 0.39 GHz, on 
average). Flux densities of the individual channels are 
plotted with 1 a error bars in Figures [2] and [3] 

While the unblocked aperture of the GBT mitigates 
the significance of side lobes, elevation scans of the cal- 
ibration source showed side lobes with an amplitude of 
«3% of the beam peak. The added power may cause 
an overestimation of the 33 GHz flux density. Using a 
model beam having side lobes at 5% of the beam peak, 
we recalculated the photometry for each radio, submil- 
limeter, and infrared image finding that our results are 
not significantly affected by the presence of side lobes 
(see §3). 

Since the reference beam throw is only 1'3, real signal 
in our reference positions is a concern; having reference 
positions on the galaxy may result in an underestima- 
tion of the 33 GHz flux density (see bottom panel of 
Figure [J). To assess the severity of this, we add back to 
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TABLE 1 

Positions, Ka-BAND Flux Densities, and Excesses 







R.A. 


Decl. 


^33 GHz 


Jexc 

•> 33 GHz 




ID 




(J2000) 


(J2000) 


(mjy) 


(%) 


corrt 


Nucleus 


20 34 52.34 


+60 9 14.2 


15.7+0.79 


-26+16 


1.00 


Enuc. 


1 


20 35 16.65 


+60 11 1.1 


0.4+0.05 


-137+47 


2.57 


Enuc. 


2 


20 35 25.49 


+60 10 1.8 


2.4+0.16 


10+11 


1.05 


Enuc. 


3 


20 34 51.89 


+60 12 44.8 


1.0+0.07 


14+15 


1.08 


Enuc. 


i 


20 34 19.17 


+60 10 8.7 


2.9+0.15 


45+ 6 


1.06 


Enuc. 


5 


20 34 39.27 


+60 4 55.1 


0.5+0.05 


-22+28 


1.19 


Enuc. 


6 


20 35 6.09 


+60 11 0.6 


2.8+0.18 


-9+13 


1.09 


Enuc. 


7 


20 35 11.21 


+60 8 59.7 


3.1+0.20 


10+10 


1.03 


Enuc. 


8 


20 34 32.52 


+60 10 22.0 


1.8+0.14 


-18+15 


1.26 


Enuc. 


9 


20 35 12.71 


+60 8 52.8 


2.4+0.16 


8+11 


1.05 



Note. — t Estimated S33 GHz correction factor due to over sub- 
traction by reference beams. These corrections were not included in 
the calculation of /3f C Q Hz - 



the 33 GHz flux densities the average flux densities mea- 
sured among the reference positions using two modeled 
33 G Hz maps by: 1) scalin g the 24 /im data (Equation 
10 in IMurphv et aI1[200cl . and 2) scaling the 8.5 GHz 
data with the measured spectral indices from fitting the 
radio data shortward of 10 GHz. The minimum flux 
density of the reference beam positions for each target 
was then taken as the local sky and subtracted. Taking 
the median value among all reference positions appearing 
clearly off the disk as a global sky measurement yields 
very similar results. These 'boosting' factors, taken as 
the weighted average between the values obtained using 
the two different modeled 33 GHz maps, are listed in 
Tableland are typically small (<10%) with the excep- 
tion of extranuclear regions 1, 5, and 8, for which one 
of the reference beams appeared to fall on an adjacent 
bright region. While the 8.5 GHz-derived 33 GHz flux 
densities at the reference positions are typically twice as 
large as the 24 /im-derived values, which only estimate 
the free-free contribution, the corresponding 'boosting' 
factors agree within 5%, on average (i.e. 1.04 ± 0.05). 

2.4. Modeling the Infrared-to- Radio Spectra 

We fit the infrared and radio data indepen dently us- 
ing th ree separate components as described in [Murphy 
(2009); thermal dust emission, thermal (free-free) ra- 
dio emission, and non-thermal synchrotron emission (see 
Figures [2] and [3]) . The infrared, and submillimeter data 
where available , wer e fit using the SED templates of 
IDale fc Heloul (|2002D . We then fit the radio data us- 
ing a combination of thermal and non-thermal emission. 
Since we are interested in measuring excess emission in 
the microwave range, we perform two separate fits using: 
1) all (1.4 — 40 GHz) radio data, and 2) only radio points 
at frequencies below 10 GHz since these frequencies are 
likely free of significant anomalous dust emission. 

We first assume a thermal (free-free) component which 
scales as S*J oc v~ aT , where q?t = 0.1 is the thermal 
spectral index. The non-thermal component is propor- 
tional to the energy loss of CR electrons to synchrotron 
emission relative to the total energy loss, scaling as 
S* ,T oc \dE/dt\ syn N(E) oc ^" aNT , where N(E) is the 
number density of CR electrons per unit energy and «nt 
is the non-thermal spectral index. We assume an en- 
ergy injection spectrum of Q(E) — nE~ p , where p = 2.8 



for shock-accelerated electro ns that have interacted with 
the ambient med ium (e.g. iBiermarm fc Strom I [19931: 
iBecker et al. 1 120091) . Our models of the non-thermal 
emission formally include energy losses to CR electrons 
arising from synchrotron, inverse Compton (IC) scat- 
tering, ionization, bremsstrahlung, as well as a pre- 
scrip tion for escape ; losses from advection are ignored 
(see IMurphv 1 [2009. for details). For simplicity, we as- 
sume interstellar medium (ISM) densities of nisM = 0.1 
and 10 cm -3 for the extranuclear regions and nucleus, re- 
spectively. Consequently, ionization and bremsstrahlung 
losses are insignificant relative to synchrotron and IC 
losses for the extranuclear regions, and become compa- 
rable for the nucleus. 

Under these assumptions we generate the shape of 
the non-thermal radio spectrum, resulting in values of 
Qnt ~ 0.8, on average for the extranuclear regions, and 
«nt ~ 0.7 for the nuclear region, and vary the fraction 
of thermal emission by scaling the relative amplitude of 
the thermal and non-thermal components to best fit all 
radio points, and only those below 10 GHz, using a stan- 
dard x 2 minimization procedure (dot dashed- and thin 
solid-lines in Figures [2] and [3l respectively) . The corre- 
sponding best-fit model 33 GHz thermal fractions, rang- 
ing between w 50 — 85%, are given in the top right corner 
of each panel along with their uncertainties. For refer- 
ence, the global thermal fraction at w 33 GHz is roughl y 
50% among star- forming galaxies (|Condon fc Yin 1 11990). 
Errors in the fitting and associated parameters are cal- 
culated by a standard Monte Carlo approach using the 
photometric uncertainties of the input flux densities. 

3. RESULTS: EXCESS EMISSION AT 33 GHZ 

We find that the Ka-band continuum is typically dom- 
inated by free-free emission, having a median thermal 
fraction of «70%. We also estimate whether there is 
33 GHz emission in excess to what our fits to the lower 
frequency (y < 10 GHz) radio data would predict (i.e. 

/33°GHz = (^33 S GHz — ^33°GHz)/^33 S ghz) an d attach an 
uncertainty (see Table 1). Half of the sample have posi- 
tive excesses ranging between w 10 — 50% (i.e. extranu- 
clear regions 2, 3, 4, 7, and 9), but only one (extranuclear 
region 4) has a statistically significant excess (« 7a). 
From our fitting, we find that «50% of the observed flux 
density from this region arises from free-free emission 
while the other fti50% appears to be anomalous, consis- 
tent with iDickinson et al. 1 (|2009f ) who report compara- 
ble fractions of free- free and anomalous dust emission for 
two components of a Galactic H11 region. Even when us- 
ing the photometry obtained with the exaggerated side- 
lobes (see §2.3), the excess in extranuclear region 4 is still 
detected at the « 5cr level indicating the detecti on is ro- 
bust. We also apply Peirce's criterion for outliers (jPeirce I 
Il852f ) to the distribution of all excesses, corrected for the 
over subtraction of sky (§2.3), finding extranuclear region 
4 to be an outlier at a significance of ~5a pop . 

For extranuclear region 1 we find that the measured 
flux density is significantly lower (i.e. «40%) than pre- 
dicted; we attribute this discrepancy to over subtraction 
of sky as one of the reference beams fell on extranuclear 
region 6. The corrected flux density, given by a cross in 
Figure [2l is very close to the expected value. 

The choice of p for the electron injection spectrum af- 
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Fig. 2. — The radio spectrum of each st ar-forming region excluding extranuclear region 4 (see Figure [3). The best-fit model spectra 
are given by thick solid lines. The Dale & Hclou (2002) templates which best fit the infrared and submillimeter data are shown using the 
short dashed-lines. The free-free, non-thermal, and spinning dust components of the radio spectra are shown using triple dot dashed-, long 
dashed-, and dotted-lines, respectively. The best-fit models neglecting the anomalous dust c omp onent which fit all the radio data, and those 
data below 10 GHz, are given by the dot dashed- and thin solid-lines, respectively (See £|2,4|l . For the nucleus and extranuclear regions 
1, 5, 6, and 8, the best-fit models overlap the dot dashed-lines since excess 33 GHz emission was not detected. In the top right corner of 
each panel the 33 GHz thermal fraction from the best-fits are given along with the best-fi t ISM densities and kinetic gas temperature for 
the spinning dust component under WNM or WIM conditions (Draine & Lazarian 1998b) where appropriate. For extranuclear regions 1, 
5, and 8, where we believe a considerable amount of signal was subtracted out by the reference beams, a cross is placed at the position 
of the corrected flux density. Below each spectrum the signal-to-noise of the residuals are shown for the two- (dot- dashed-lines) and 
three-component (thick-solid-lines) models along with corresponding x 2 values (see §4.1). 
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Fig. 3. — The radio spectrum of extranuclear region 4 along with 
signal-to-noise of the residuals for the best-fit models (see caption 
of Figure In- 



fects the average 33 GHz thermal fraction from the fits, 
along with the significance of the excesses. For exam- 
ple, choosing p w 2.2, which is in the injection limit 
(2.0 < p < 2.4) for electrons freshly shock-accelerated 
by SN remnants, yields an average thermal fraction at 
33 GHz of <10% and results in the excesses from regions 
2, 3, 7 and 9 also being significant at >2.5cx level. How- 
ever, the x 2 values from these fits are >3 times larger, on 
average, suggesting that an injection-limit value for p is 
inappropriate. On the other hand, setting p ks 2.6 gives 
similar \ 2 values and a median 33 GHz thermal fraction 
of w60%. 

4. DISCUSSION: THE NATURE OF THE 33 GHZ 
EXCESS 

We report on the likely detection of excess microwave 
emission from an extranuclear star-forming region in 
NGC 6946. A rising spectrum at these frequencies can be 
produced by optically thick free-free emission observed 
from compact Hn regions. However, this explanation 
appears inconsistent given the constraints set on the pro- 
duction rate of Lyman continuum photons (AL yc ) by the 
24 //m data. Assuming that all of the 33 GHz flux density 
from extranuclear region 4 arises from free-free emission 
suggests Al Y c ~ 1.4 x 10 52 s _1 , a factor of more than 
5 times larger than the N-^ yc ~ 2.6 x 10 51 s _1 rate in- 
ferred from the 24 /Ltm data. Thus, explaining the excess 
33 GHz emission by compact Hn regions seems unlikely. 

It is at precisely these frequencies where excess emis- 
sion, coined 'anomalous' dust emiss ion, has been found 
in diffuse emission e xperiments fe.g.lDavies et al. Il2006t 
iHildebrandt et al. I l2007t iDobler fc Finkbeiner I l2008f ). 
Taking extranuclear region 4 as an example, we find 
that its 33 GHz flux density is better fit by including 
a component besides free-free, synchrotron and thermal 
dust emission, and its shape appears reasonably well con- 
strained (see Figure [3]) . On the high frequency side, the 
spectrum is constrained by measurements from the in- 
frared into the submillimeter, and is most unlikely to 



exceed the iDale fc Helou I (|2002[ ) model curves above 
300 GHz. On the low frequency side, data from extranu- 
clear region 4 requires the excess emission component to 
fall off from 30 to 10 GHz. This supports the spinning 
dust interpretation, which predicts a peaked component 
in this frequency range. We therefore believe that our 
observations indicate that an anomalous dust component 
exists, it is detectable in extragalactic sources, and it can 
be reproduced by spinning dust. To our knowledge, this 
is the first detection of anomalous dust emission outside 
the Galaxy. 

4.1. The Spinning Dust Model 

Currently, the favored origin of this excess emission is 
electric dipole radiation from rapid ly rotating ultrasmall 
grains dDraine fc Lazarian I Il998bf). Using t he rev ised 
model of A li-Haimoud. Hirata. fc Dickinson I (|2009f ). we 
fit all occurrences of excess Ka-band emission (see §3) as- 
suming a range of warm neutral medium (WNM; sh ~ 
0.01) and warm ionized medium (WIM; x-r sa 0.99 ) 
conditions as defined by IDraine fc Lazarianl ([l998b), 
where xn is the ionization fraction. The parameters 
that were allowed to vary were ISM density (nisM = 
0.1, 0.3, 0.5, 1, 5, 10, 15, 30 cm~ 3 ) and gas kinetic temper- 
ature (T = 0.1,0.3,0.5,0.7,1 x 10 4 K). The spinning 
dust spectra were first scaled to match the observed ex- 
cess 33 GHz flux density and then added to the best- 
fit thermal+non-thermal radio spectra from fitting the 
radio data below 10 GHz. The residuals of this three 
component model against the data, fitting each of the 4 
CCB frequency channels independently, were minimized 
to pick the favored spinning dust component by varying 
the spinning dust parameters only. The signal-to-noise 
of the residuals, along with the corresponding x 2 values, 
are shown beneath each radio spectrum in Figure for 
the cases of fitting all radio data with two components 
(i.e. thermal + non-thermal radio emission) and this 
three component model. 

The emissivity of the best-fit spinning dust model for 
extranuclear region 4 corresponds to a hydrogen column 
density of Nh ~ 3.3 x 10 21 cm -2 . This is a factor of 
«1.7 times larger than the Nh ~ 1.9 x 10 21 cm" 2 aver- 
age co lumn measured using Hi data from IWalter et al~l 
(2008), and therefore consistent with the observed levels 
for a reasonable molecular gas fraction. 

While the excess at 33 GHz is only statistically signif- 
icant for extranuclear region 4, the fits are improved by 
including the spinning dust component in extranuclear 
regions 2, 3, 7, and 9 as evidenced by the lower x 2 val- 
ues. The ISM density and gas temperature from the best 
fitting spinning dust component are given and point to 
WNM parameters for each region. However, the lack of a 
fine spectral sampling, and large uncertainty in the Ka- 
band spectral index, yields a large range of acceptable 
values. It should also be noted that each of these regions 
are located near the disk-edge of NGC 6946. 

4.2. Implications 

These results suggest that the Ka-band observations 
from external galaxies may contain appreciable amounts 
of anomalous dust emission and that the assumption that 
this frequency range is dominated by free-free emission, 
making it a clean measure for the star formation in exter- 
nal galaxies, needs further investigation. It is currently 
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unclear how much of an effect this component might have 
for interpreting globally integrated 33 GHz flux density 
measurements of galaxies. Summing the 33 GHz flux 
densities for all 10 regions, we find that the excess emis- 
sion accounts for <10% of the total 33 GHz flux density 
even when correcting for badly placed reference positions. 
Future Ka-band observations at higher resolution using 
the EVLA should help to address this by precisely iden- 
tifying the origin of the excess microwave emission (e.g. 
in the PDRs surrounding the Hn regions), and in turn 
also allow for signifiant progress in modeling this com- 
ponent out of the Galactic foreground emission for CMB 
experiments. 
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